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We use a combination of numerical and analytical methods, exploiting the equations derived in a 
preceding paper, to classify all spherically symmetric self-similar solutions which are asymptotically 
Friedmann at large distances and contain a perfect fluid with equation of state p = (7 — l)fJ. with 
< 7 < 2/3. The expansion of the Friedmann universe is accelerated in this case. We find a 
one-parameter family of self-similar solutions representing a black hole embedded in a Friedmann 
background. This suggests that, in contrast to the positive pressure case, black holes in a universe 
with dark energy can grow as fast as the Hubble horizon if they are not too large. There are also 
self-similar solutions which contain a central naked singularity with negative mass and solutions 
which represent a Friedmann universe connected to either another Friedmann universe or some 
other cosmological model. The latter are interpreted as self-similar cosmological white hole or 
wormhole solutions. The throats of these wormholes are defined as two-dimensional spheres with 
minimal area on a spacelike hypersurface and they are all non-traversable because of the absence of 
a past null infinity. 

PACS numbers: 04.70.Bw, 95.36.+X, 97.60.Lf, 04.40.Nr, 04.25. Dm 



I. INTRODUCTION 

The study of black holes in stationary and asymptoti- 
cally flat spacetimes has led to many remarkable insights, 
such as the "no hair" theorem and the discovery of the 
thermodynamic properties of black holes. In contrast, 
black holes in a Friedmann background have not been 
fully investigated and many important open questions 
remain. When the size of the black hole is much smaller 
than the cosmological horizon, asymptotically flat black 
hole solutions will presumably provide a good approxima- 
tion. However, for black holes comparable to the horizon 
size - such as primordial black holes (PBHs) formed in 
the early universe [l| - one would need to take into ac- 
count the effects of the cosmological expansion. 

The main issue addressed in this paper is the growth 
of such "cosmological" black holes. Zel'dovich and 
Novikov flrst discussed the possibility of self-similar 
growth for a black hole in a Friedmann universe in 1967. 
This suggested that the size of a PBH might always be 
the same fraction of the size of the particle horizon until 
the end of the radiation-dominated era, when they would 
have a mass of order IO^^Mq. Since there is no evidence 
that such huge black holes exist, this seemed to imply 
that PBHs never formed at all. More recently, the possi- 
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bility of self-similar PBH growth in the quintessence sce- 
nario has led to the suggestion that PBHs could provide 
seeds for the supermassive black holes thought to reside 
in galactic nuclei However, both these arguments are 
based on a simple Newtonian analysis. 

The first general relativistic analyses of self-similar cos- 
mological black holes was given by Carr and Hawking [3| 
for a radiation fluid in 1974. They showed that there is 
no self-similar solution which contains a black hole at- 
tached to an exact Friedmann background via a sonic 
point (i.e. in which the black hole forms by purely causal 
processes). This result was extended to the perfect fluid 
case with p = (7 — l)/i and 1 < 7 < 2 by Carr Q and 
Bicknell and Henriksen Q. For a stiff fluid (7 = 2), it 
was originally claimed by Lin et al. 0] that self-similar 
cosmological black holes could exist. However, Bicknell 
and Henriksen 0] subsequently showed that this solution 
does not represent a black hole after all, since Lin et al. 
had misidentificd the event horizon. Although Bicknell 
and Henriksen did construct some numerical self-similar 
solutions, in these the stiff fluid turns into null dust at a 
timelike hypersurface, which seems physically implausi- 
ble. Indeed, we have recently proved rigorously the non- 
existence of self-similar cosmological black holes for both 
a stiff fluid and a scalar field Q . 

In the present paper, we focus on fluids with < 7 < 
2/3. Although this equation of state violates the strong 
energy conditions, it may still be very relevant for cos- 
mology - both in the early universe (when inflation oc- 
curred [l3l) and at the current epoch (when the acceler- 
ation may be driven by some form of dark energy pd|). 
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Such matter exhibits "anti-gravity" - in the sense that 
the active gravitational mass in the Raychaudhuri equa- 
tion is negative - but the dominant, null and weak energy 
conditions still hold. Recently, the accretion of dark en- 
ergy or aphantom field onto a black hole has also been 
studied ITsj l , but the cosmological expansion is again 
neglected in these analyses. 

The < 7 < 2/3 case is particularly important in 
the cosmological context because there then exists a one- 
parameter family of self-similar solutions which are ex- 
actly asynrptotic to the flat Friedmann model at large 
distances In the positive pressure case, the solu- 

tions cannot be "properly" asymptotic Friedmann be- 
cause they exhibit a solid angle deficit at infinity, so 
it would be more accurate to describe them by "quasi- 
Friedmann" p^ . In the present paper, we obtain a one- 
parameter family of physically reasonable self-similar cos- 
mological black hole solutions numerically. This strongly 
suggests that a black hole in a universe filled with dark 
energy can grow in a self-similar manner. In ref. this 
system was also numerically investigated but no black 
hole solution was reported. This is because one needs to 
consider the analytic extension of the solutions to find 
the black holes and this was not done in ref. [l^l • 

We also find that there exists a class of self-similar 
cosmological wormhole solutions. A wormhole is an 
object connecting two (or more) infinities. Einstein's 
equations certainly permit such solutions. For exam- 
ple, the well-known static wormhole studied by Morris 
and Thorne [1^1 connects two asymptotically flat space- 
times. Although the concept of a wormhole is originally 
topological and global, it is possible to define it locally 
by a two-dimensional surface of minimal area on a non- 
timelike hypersurface. However, in any local definition, 
we inevitably face the problem of time-slicing, because 
the concept of "minimal area" is slice-dependent. 

Hochberg and Visser and Hayward [ll] define a 
wormhole throat on a null hypersurface, in which case, 
a light ray can travel from one infinity to another infin- 
ity. Their definition excludes the maximally extended 
Schwarzschild spacetime from the family of wormhole 
spacetimes. It also implies that in the non-static situ- 
ation the null energy condition (/.( > for the present 
matter model) must be violated at the wormhole throat, 
so the existence of wormholcs might seem implausi- 
ble [13, S, [l^- However, their definition of a wormhole 
is physically reasonable only in the presence of a past 
null infinity, and in the cosmological situation this often 
does not apply because there is a big-bang singularity. 
Thus, wormholes in their sense may be too restrictive in 
the cosmological context. 

This motivates us to define a wormhole throat on a 
spacelike hypersurface. Here we adopt the conventional 
definition of traversability in which an observer can go 
from one infinity to another distinct infinity. Our numer- 
ical solutions are not wormholes in the sense of Hochberg 
and Visser or Hayward [l^ but they are wormholes 
in the sense that they connect two infinities. Also, al- 



though dark energy violates the strong energy condition, 
the null energy condition is still satisfied in our worm- 
hole solutions. Numerical solutions containing a super- 
horizon-size black hole in a Friedmann universe filled with 
a massless scalar field provide another interesting exam- 
ple of this sort of wormhole [20|. These involve a worm- 
hole structure with two distinct null infinities, where all 
possible energy conditions are satisfied. 

We find that there are three types of cosmological 
self-similar wormhole solutions. The first connects two 
Friedmann universes, the second a Friedmann and quasi- 
Friedmann universe, the third a Friedmann universe and 
a quasi-static spacelike infinity. In fact, solutions of 
the third type are naturally interpreted as a Friedmann 
universe emergent from a white hole. In order to ob- 
tain these results, we utilize the formulation and asymp- 
totic analyses presented in the accompanying paper |2l| 
(henceforth Paper I). The combination of these two pa- 
pers completes the classification of asymptotically Fried- 
mann spherically symmetric self-similar solutions with 
dark energy 

The plan of this paper is as follows. In Section |TT1 we 
present the basic variables and describe the numerical 
results. (Paper I provides the detailed formulation and 
field equations.) In Section lllll we classify the self-similar 
solutions in terms of their physical properties. Section |TV] 
provides a summary and discussion. 



II. NUMERICAL RESULTS 

We consider a spherically symmetric self-similar space- 
time with a perfect fluid and adopt comoving coordinates. 
The line element can be written as 



(2.1) 



where z = r/t is the self-similar variable, R = rS is the 
arcal radius, and we adopt units with c = 1. The Einstein 
equations imply that the pressure p, energy density fi and 
the Misner-Sharp mass m must have the form 



2Gm 



W{z) 
rM{z). 



(2.2) 

(2.3) 
(2.4) 



We assume an equation of state p = (7 — l)/.i with 
0<7<2/3. A crucial role is played by the velocity 
function V = |2;|e*~*, which is the velocity of the fluid 
relative to the similarity surface of constant z. A simi- 
larity surface is spacelike for (1 — F^)e^* < 0, timelike 
for (1 - y2)e2* > 0, and nuU for (1 - V^)e^'^ = 0. In 
the last case, the surface is called a "similarity horizon" . 
We also use trapping horizon notation for describing the 
physical properties of solutions. (See refs. [13, HI] for the 
definition of a trapping horizon.) 
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The basic field equations are derived from the Einstein 
equations. Because of the self-similarity, they reduce to a 
set of three ordinary differential equations (ODEs) for S, 
S' and W, together with a constraint equation. (This is 
discussed in Paper I and also ref. [gj] ) . Another formula- 
tion of the field equations, involving a set of three ODEs 
for M, S and W together with a constraint equation, is 
given in ref. [25| . 

We numerically evolve the ODEs given by Eqs. (2.27) 
and (2.28) in Paper I (henceforth referred to as (1.2.27) 
and (1.2.28)) from small positive z (which corresponds to 
large physical distance). We use asymptotically Fried- 
mann solutions in the form given by Eqs. (1.4.1) and 
(1.4.2), where Eqs. (1.4.13) and (1.4.14) specify the form 
of the perturbations. We set the gauge constants qq and 
bo in Eqs. (1.3.5) and (1.3.6) to be 1 for the Friedmann so- 
lution and this determines the gauge constants cq and ci 
through Eqs. (1.3.11) and (1.3.12). A single free parame- 
ter Aq then characterizes these solutions. Aq represents 
the deviation of the energy density from that of the Fried- 
mann universe at spatial infinity, with positive (negative) 
Aq corresponding to an overdensity (underdensity) there. 
We assume 7 = 1/3 as a typical example and numerically 
integrate the ODEs using the constraint equation (1.2.31) 
or (1.2.32) to monitor the numerical accuracy. 

As shown in Paper I, some of the solutions admit an 
extension beyond z = 00 (which corresponds to a finite 
non-zero physical distance). To assure this, we choose 
suitable dependent variables and then expand these as 
Taylor series with respect to a new independent vari- 
able around z = ±00. Paper I specifies the variables 
which satisfy this requirement for asymptotically quasi- 
Kantowski-Sachs and quasi-static solutions. This is our 
key technique for obtaining the extension of the solu- 
tion while retaining good numerical accuracy. In order 
to display the solutions, we plot — 1/z on the horizon- 
tal axis. The slightly different set of equations used by 
Harada and Maeda [2^ are equivalent to those given in 
Paper I. In fact, we have constructed numerical codes 
to evolve both sets of ODEs independently and find ex- 
cellent agreement. Our numerical accuracy is checked 
by the constraint equation, which is always satisfied to 
within a factor of 10~^. 

We show the numerical solutions in Figs. [T]^ Al- 
though we have used many other values for the parame- 
ter Aq, we only show the results for Aq = —0.1, —0.08, 
-0.06, -0.0404, -0.03, -0.02, -0.01, 0, 0.01 and 0.02. 
Aq — corresponds to the exact Friedmann solution. 

Figure [1] shows the velocity function V. The thick line 
corresponds to the exact Friedmann solution and this 
has one similarity horizon. The solutions with positive 
Aq also have one similarity horizon (where V = I) and 
encounter curvature singularities with negative mass at 
finite positive z, with V converging to 1 from below. The 
solutions with —0.0253 < Aq < have two similarity 
horizons and encounter curvature singularities with posi- 
tive mass at finite negative z , with V converging to 1 from 
above. The solutions with -0.0780 < Aq < -0.0253 also 
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FIG. 1: V is plotted against for 7 = 1/3. The thick line 

corresponds to the exact Friedmann solution (ylo = 0) , which 
has one similarity horizon. 



have two similarity horizons. However, of these only the 
solution with Aq 2± —0.0404 converges to the asymptoti- 
cally Friedmann solution as z — > 0~; the other solutions 
are asymptotically quasi-Friedmann. The solutions with 
Aq < —0.0780 encounter curvature singularities with pos- 
itive mass at finite negative z, with V converging to 
1 from above. The number of similarity horizons de- 
pends on the value of Aq: there are zero, one and two 
similarity horizons for Aq < —0.108, ^0 — —0.108 and 
-0.108 <Aq< -0.0780, respectively. 

Figure [2] shows the functions W = SnGfir^ and 
W/z'^ = SnGfit'^. The thick line W/z'^ = 12 corresponds 
to the exact Friedmann solution {Aq = 0) and has a 
big-bang singularity as z — > -|-oo and t ^ 0. W j z^ 
diverges at finite positive z for solutions with positive 
Aq, corresponding to the negative-mass curvature sin- 
gularities seen above. It also diverges at finite negative 
z for solutions with —0.0253 < ^0 < 0, corresponding 
to the positive-mass singularities. The solutions with 
-0.0780 < Aq< -0.0253 have no singularity in l^/z^; 
instead they converge to 12 as z — > 0~, corresponding to 
the asymptotically quasi-Friedmann solutions. 

Figure [3] shows the functions S = R/r and zS. Note 
that decreasing (increasing) S" as a function of — 1/z 
means that the spacetime is expanding (collapsing) . The 
thick line z'^S = 1 corresponds to the exact Friedmann 
solution. The curvature singularity at finite non-zero z 
in the solution with Aq > —0.0253 is a central singular- 
ity. There exists a collapsing region near this curvature 
singularity in the solution with —0.0253 < Aq < 0. By 
contrast, there exists a single local minimum in the pro- 
file of zS for the solution with Aq < —0.0253 and this 
corresponds to a wormhole throat. The solution with 
-0.0780 < Ao < -0.0253 has two different spacelike in- 
finities at z = O"*" and 0~, where zS diverges to -|-oo. The 
solution with ^0 ^ —0.0780 also has two different space- 
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FIG. 3: (a) S = R/r and (b) zS = R/t are plotted for 
FIG. 2: (a) W = SnGur'^ and (b) W/z^ = SnGfit'^ are plot- 7 = 1/3. The thick line z'^S = 1 corresponds to the exact 
ted against —1/z for 7 = 1/3. The thick line W/z'^ = 12 Friedmann solution {Ao — 0). 
corresponds to the exact Friedmann solution {Ao — 0). 



like infinities at z = 0+ and z = +00. In the latter class, 
zS diverges but S is finite, so the solution extends into 
the negative z (negative t) region, where S is negative 
for positive R. 

Figure [4] shows the functions M ~ 2Gm/r and zM — 
2Gm/t. The thick line z'^M — 4 corresponds to the 
exact Friedmann solution. We can see that the curva- 
ture singularity at finite non-zero z in the solutions with 
-0.0253 < ylo < or < -0.0780 has positive mass, 
while that in the solution with positive Aq has negative 
mass. Note that zM is negative for positive m in the so- 
lutions with < —0.0780, since negative z corresponds 
to negative t. 

Figure El shows the function M/S = 2Gm/R. Note 
that M/S > 1 and M/S < 1 correspond to the trapped 
and untrapped regions, respectively, and M/S = 1 in 



the collapsing (expanding) region corresponds to a fu- 
ture (past) trapping horizon. The thick line corresponds 
to the exact Friedmann solution; this has a single past 
trapping horizon, corresponding to a cosmological trap- 
ping horizon. The solutions with positive Aq also have a 
single past trapping horizon. Those with Aq < —0.0253 
have no trapping horizon and the spacetime is trapped 
everywhere. Those with —0.0253 < < have one past 
trapping horizon and one future trapping horizon. The 
latter corresponds to a black hole trapping horizon. 

From these numerical results, we classify asymptoti- 
cally Friedmann solutions as z — > 0^ into six classes: 
(a) the flat Friedmann universe (^o = 0); (b) a cos- 
mological naked singularity {Aq > 0); (c) a cosmologi- 
cal black hole (ai < Aq < 0); (d) a Friedmann-quasi- 
Fricdmann wormhole (as < Aq < ai with Aq ^ a2); (e) 
a Friedmann- Friedmann wormhole {Aq = 012); and (f) 
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FIG. 4: (a) M = 2Gm/r and (b) zM = 2Gm/t are plotted 
for 7 = 1/3. The thick hne z^M = 4 corresponds to the exact 
Friedmann solution (j4o = 0). 



a cosniological white hole (/lo < as)- For the present 
choice of gauge constants, a\ ~ —0.0253, ol2 — —0.0404 
and as - -0.0780. 

The solutions in class (b) describe a negative mass 
naked singularity in a Friedmann universe. Those in class 
(c) describe a positive mass black hole in a Friedmann 
universe. Those in class (d) contain a wormhole throat 
connecting a Friedmann universe and a quasi-Friedmann 
universe. Class (e) has only one solution and this con- 
tains a wormhole throat connecting two Friedmann uni- 
verses. The solutions in class (f) describe a positive mass 
white hole in a Friedmann universe which contains a 
wormhole throat connecting a Friedmann universe and 
a quasi-static spacelikc infinity. 

The exact flat Friedmann solution, the only solution of 
class (a), is located at the threshold between classes (b) 
and (c), while we could not resolve the threshold solution 




FIG. 5; M/S = 2Gm/R is plotted for 7 = 1/3. The thick 

0.02 line corresponds to the exact Friedmann solution (Ao — 0). 

0.01 



-0.01 

-0.02 

).0404 between classes (d) and (f). The solution of class (c) is 

008 ------- obtained by fine-tuning the parameter A^. We investi- 

-0.1 gate their properties and physically interpret them in the 

next section. 



III. PHYSICAL INTERPRETATION 



A. Exact flat Friedmann universe 



z — -|-0, t = 00 




FIG. 6: The conformal diagram of the exact flat Friedmann 
solution. The thin curves and lines denote similarity surfaces, 
i.e., orbits of z ^constant. There is a similarity horizon at 
z = zi (> 0). 



The solution with = is the exact flat Friedmann 
model. The limit z +00 corresponds to the physical 
center i? = at constant t or to the big-bang singularity 
at constant r. The conformal diagram of this spacetime 
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A Friedmann 
^ infinity 
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FIG. 7: The t-r plane of the exact flat Friedmann solution. 
A fluid element moves from the bottom to the top. r = and 
r — > +00 correspond to spacelike infinity and the physical 
center (7? = 0), respectively, t = corresponds to the big- 
bang singularity, which is represented by a zig-zag line. There 
is a similarity horizon at z = z\ (> 0), corresponding to a 
cosmological event horizon. 



Friedmann 
^/infinity 




FIG. 9: The t-r plane of the cosmological naked singular- 
ity solution, r = corresponds to spacelike infinity. There 
is a similarity horizon eX z = z\ (> 0), corresponding to a 
cosmological event horizon, z = z* (> 0) with < i < cxd 
corresponds to a timelike singularity of negative mass and is 
represented by a zigzag line. 



and the t-r plane of the solution are shown in Figs. [5] and 
[71 respectively. There is a similarity horizon at z = zi, 
where V = 1, corresponding to the cosmological event 
horizon and the initial singularity is null. 

B. Cosmological naked singularity 



z = +0, t = oo 



From Fig. [31 we see that R,r > everywhere and also 
R^t > (i.e. the spacctime is expanding), where a comma 
denotes a partial derivative. Figures [1] and [5] show that 
there is a cosmological event horizon and a cosmological 
trapping horizon. Near the singularity, V approaches 1 
from below according to Eq. (1.4.60) with negative Vi. 
Figures [4] and \5\ show that the singularity has negative 
mass and is in the untrapped region. This solution de- 
scribes a naked singularity embedded in a Friedmann 
background. The conformal diagram of the spacctime 
and the t-r plane of the solution are shown in Figs. [8] and 
[HI respectively. There is a similarity horizon at z = zi 
and this corresponds to the cosmological event horizon. 



FIG. 8: The conformal diagram of the cosmological naked 
singularity solution. There is a similarity horizon at z = zi 
(> 0). z = 2* (> 0) with < t < oo corresponds to a timelike 
singularity at which the mass is negative. 



The solutions with Aq > represent asymptotically 
negative-mass singularities as z ^ z* (0 < z* < oo). 



C. Cosmological black hole 

The solutions with ai < Aq < (ai ~ -0.0253) 
are asymptotically quasi-Kantowski-Sachs as z — ^ ±oo 
and asymptotically positive-mass singular as z — > z* 
(— oo < z* < 0). Because the solutions are asymptoti- 
cally quasi-Kantowski-Sachs as z — > -f oo, they can be ex- 
tended analytically into the region with negative z. Since 
the Kantowski-Sachs solution has a curvature singularity 
at f = 0, the analytic extension in this case must be in- 
terpreted as an extension from the positive r to negative 
r region. In the extended region, S and AI are negative 
but the areal radius R = rS and mass m = rM / (2G) are 
positive because r is negative. 

From Fig. [31 we find that R^r > holds everywhere 
and there exists a collapsing region (with R t < 0) near 
the singularity. It is seen in Fig. [5l that there are two 
trapping horizons, one in the expanding region and the 
other in the collapsing region. These arc the cosmological 
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and black hole trapping horizons, respectively. Figure [TOl 
shows that the trapping horizon is degenerate for the crit- 
ical value Aq = Q!i. Figures [D or [TOl show that there is 
both a cosmological and black hole event horizons. For 
Aq = —0.01 and —0.02, we see that the black hole and 
cosmological event horizons are in the untrapped and 
trapped regions, respectively. Near the singularity, V 
approaches 1 from above according to Eq. (1.4.60) with 
positive Vi . Figures [4] and [5] show that the singularity 
has positive mass and is in the trapped region. 

The forms of zS = R/t for the various horizons are 
plotted as functions of the parameter in Fig- [HI The 
ratio of the black hole event horizon size to the cosmo- 
logical event horizon size goes from to 0.36, while the 
ratio of the black hole event horizon size to the back- 
ground Hubble horizon size goes from to 0.70. On the 
other hand, the ratio of the cosmological trapping hori- 
zon size to the black hole trapping horizon size goes from 
to roughly 1. 

This solution describes a black hole in a Friedmann 
background. Figures fT2l and [T3l show the conformal dia- 
gram and t-r plane for the cosmological black hole solu- 
tion. We can see that the initial singularity is null and 
that there exists both a null and spacelike portion of the 
black hole singularity. There are two similarity horizons 
at zi > and Z2 < 0. These correspond to the cosmo- 
logical and black hole event horizons, respectively. 



D. Friedmann-quasi- Friedmann wormhole 

The solutions with < Aq < a\ (03 ~ —0.0780) are 
also asymptotically quasi-Kantowski-Sachs as z — > ±00 
and the form of these solutions in the extended region 
can again be obtained numerically. Figure [T3] shows that 



12 and z^S ^ constant as 



0' 



In the 

asymptotically Friedmann case (Aq — 012 — —0.0404) 
the constant is —1. Otherwise the solutions are asymp- 
totically quasi-Friedmann as z — > 0~ and the constant is 
different from —1. 

From the profile of zS in Fig. [3l we see that a solution 
in this class has a single throat at z = zt, characterized 
by R^r = and R,r < (>)0 for -1/z < (>) - 1/zt, i.e. 
the areal radius R has a minimum on the spacelike hy- 
persurface with constant t. As a result, the spacetime 
represents a dynamical wormhole with a throat connect- 
ing a Friedmann universe to another Friedmann or quasi- 
Friedmann universe. With Aq = a2, the solution has re- 
flection symmetry, so the wormhole throat connects two 
exact Friedmann universes and is located at z = ±00. 
Fine-tuning of Aq is required because the asymptotically 
Friedmann solution given by Eqs. (1.4.13) and (1.4.14) is 
non-generic among the solutions of the linearized equa- 
tions for A, A', B' and B" in the region where V"^ ^ 00 
(see Paper I). 

This wormhole universe is born from the initial singu- 
larity at i = 0, around which the spacetime for fixed r 
is described by the quasi-Kantowski-Sachs solution. The 
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FIG. 10; (a) V and (b) M/S as functions of zS{= R/t) for 
the cosmological black hole solutions. 



profile of S in Fig. [3] shows that the spacetime is ex- 
panding everywhere and Fig. 3] shows that it has positive 
mass. From Figs. [T] and O there are two similarity hori- 
zons, corresponding to two cosmological event horizons, 
but there is no trapping horizon. 

Figures [T51 and [TBI show the conformal diagram and t-r 
plane for the Friedmann-quasi-Fricdmann wormhole so- 
lution. We can sec that the initial singularity is null and 
that there are two parts of future null infinity, both be- 
ing spacelike and having Friedmann and quasi-Friedmann 
asymptotics. There are two similarity horizons at zi > 
and Z2 < 0, both corresponding to cosmological event 
horizons. Figure [5] shows that the spacetime is trapped 
everywhere, so there is no trapping horizon. This explicit 
example shows that the wormhole definitions of Hochbcrg 
and Visser 17 1 or Hayward [l^ are too restrictive. 
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Cosmological event horizon 
Black hole event horizon 
Cosmological trapping horizon 
Black hole trapping horizon 



-0.03 -0.025 -0.02 -0.015 -0.01 -0.005 
Ao 

FIG. 11: The value of zS = R/t for various horizons is plotted 
as a function of the parameter Aq for the cosmological black 
hole solutions with 7 = 1/3. 



z = +0, t = 00 



v. 




FIG. 12: The conformal diagram of the cosmological black 
hole solution. There are two similarity horizons at z = 21 
and z = Z2 (22 < < zi). z = zi with < t < 00 is the 
cosmological event horizon, while z = Z2 with < t < 00 is 
the black hole event horizon. 2 = 2, (< 0) with < t < 00 
gives a spacelike singularity with positive mass. 
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0.01 


-0.01 
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-0.03 
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-0.06 
-0.08 
-0.1 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

-1/z 

(b) 

FIG. 14: (a) W/z^ and (b) z^S are plotted for 7 = 1/3. 



E. Cosmological white hole 




FIG. 13: The t-r plane of the cosmological black hole solution, 
r = 0"*" corresponds to spacelike infinity. 



The solutions with < are asymptotically quasi- 
static as 2; — > ±00, so that S = R/r converges to a con- 
stant there (see Fig. [3]) . We can therefore extend them 
analytically into the negative z region, which corresponds 
to negative t. 

The variables S and W can be used to obtain the solu- 
tions in the extended region because they are finite in the 
asymptotically quasi-static solutions at z = ±00. In the 
extended region, zS and zM are negative but the areal 
radius R = rS and mass m = rM/{2G) are positive be- 
cause t is negative. 

The profile of zS = R/t in Fig. [3] shows that its min- 
imum value corresponds to a throat connecting a fiat 
Fricdmann infinity and a quasi-static infinity, so this de- 
scribes a dynamical wormhole. The spacetimc is expand- 
ing everywhere and has positive mass from Fig. 31 Fig- 
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z = —0, t = oo z = +0, t = oo 




FIG. 15: The conformal diagram of the Friedmann-quasi- 
Friedmann wormhole solution. There are two similarity hori- 
zons, z = zi and z = Z2 {z2 < < zi), both corresponding to 
cosmological event horizons, z = 0^ and z — 0^ give two 
distinct null infinities, described by Friedmann and quasi- 
Friedmann asymptotics, respectively, z — ±oo, however, is 
described by the Kantowski-Sachs asymptotic as t increases 
from to oo. 




FIG. 16: The t-r plane of the Friedmann-quasi-Friedmann 
wormhole solution. We show the solution with a wormhole 
throat at 2 = zt in the positive z region, r — O'^ and 0^ 
correspond to the Friedmann and quasi-Friedmann infinities, 
respectively, t = corresponds to the big-bang singularity, 
represented by the quasi-Kantowski-Sachs solution. 



ure[5]shows that it is trapped everywhere, so there is no 
trapping horizon. 

As seen in Figs. [2] and [H there is a central curvature 
singularity at z = in the negative z region and the 
mass is positive there from Fig. This corresponds to 
the singular positive-mass asymptote discussed in Paper 
I. Near the singularity, V is approximated by Eq. (1.4.60) 
with positive Vi and approaches 1 from Fig. [TJ The sin- 
gularity is spacelike because of the positive mass and re- 
sembles a Schwarzschild white hole singularity (although 
there is no trapping horizon). This singularity is non- 
simultaneous, so it differs from the big-bang singularity. 
On the other hand, the singularity at t = with z ~ zi 
is massless and does correspond to the big-bang singular- 
ity. The solution describes a white hole in the Friedmann 



universe. 

As seen in Fig. [TJ the solutions can be classified into 
three types, according to whether they have no similar- 
ity horizon (Aq < a^), one degenerate similarity hori- 
zon {Aq = a^) OT two non-degenerate similarity horizons 
(q!4 < Aq < as), where ~ —0.108. The degenerate 
case is obtained by fine-timing Aq. 

Figures [171 and [T51 show the conformal diagram and t- 
r plane for these three types of Friedmann-quasi-static 
wormhole solutions. In the case with similarity hori- 
zon(s), the singularity has both spacelike and null por- 
tions, which correspond to the white hole and initial sin- 
gularities, respectively, and there are two parts of future 
null infinity. One is spacelike and described by the Fried- 
mann asymptote and the other is null. In the case with 
no similarity horizon, the singularity and future null in- 
finity are both spacelike. 

In closing this section, a comment should be made 
about what is meant by the term "wormhole" . It is 
sometimes assumed that wormholes require negative en- 
ergy density and there are certainly theorems to this ef- 
fect [13, [3 [ill . Our solutions do not require negative 
energy density but this does not conflict with these theo- 
rems because we have a different definition of a wormhole 
throat. 

The theorems claiming the violation of the null energy 
condition at a wormhole throat [l3, [l3l define the throat 
on a null hypersurface. In this case, the throat coincides 
with a trapping horizon. In the present paper we define 
a wormhole throat on a spacelike hypersurface. In fact, 
our wormhole solutions do not have a trapping horizon at 
all, so they are not wormholes in the sense of Hochbcrg 
and Visser [l3l or Hayward [3 even though they possess 
two distinct infinities. We discuss this problem in a more 
general framework in a separate paper [2^. 

Here we also consider the threshold between the 
Friedmann-quasi-Friedmann solutions and the cosmolog- 
ical white hole solutions. Strictly speaking, the threshold 
solution with the exact parameter value Aq = cannot 
be achieved numerically. Among the asymptotic solu- 
tions obtained in Paper I, the constant-velocity asymp- 
tote will be the most plausible candidate for this thresh- 
old solution. The resulting conformal diagram would be 
given by Fig. 1191 admitting a null ray traveling from the 
constant velocity infinity (z — > +oo,r — > +oo) to the 
Friedmann infinity (z — > +0,t — > oo). However, this pos- 
sibility is clearly excluded because this diagram contains 
a traversable dynamical wormhole, which is prohibited 
by the discussion in Hochberg and Visser [l3| or Hay- 
ward [Hi. 

IV. SUMMARY AND DISCUSSION 

We have numerically investigated spherically symmet- 
ric self-similar solutions for a perfect fluid with p = 
(7 — l)fi (0 < 7 < 2/3) which arc properly asymptotic to 
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z = +0, t = oc 




z = ±oo.t = ±0 



2 = 2» 

(b) 

z = +0, f = oo 




2 = ±0C, t = ±0 



2 = 2* 

(c) 

FIG. 17: The conformal diagrams of the cosmological white 
hole solutions with: (a) no similarity horizon; (b) one degen- 
erate similarity horizon z = z\ (> 0); and (c) two distinct 
similarity horizons z = z\ and z = Z2 {Q < Z\ < Z2 < oo). For 
each case, z = (< 0) is a spacelike singularity with positive 
mass. In (b), z = z\ with t = oo corresponds to future null 
infinity. In (c), z — Z2 with t — oo corresponds to future null 
infinity. 



the flat Friedmann model at large distances, in the sense 
that there is no solid angle deficit. We have integrated 
the field equations for the self-similar solutions and in- 
terpreted the solutions using the asymptotic analysis of 
Paper I. The results arc summarized in Tabic I. 

We have seen that there is a class of asymptotically 
Friedmann self-similar black hole solutions and this con- 
trasts to the situation with positive pressure fluids or 
scalar fields, where there are no such solutions. This 
suggests that self-similar growth of PBHs is possible in 
an inflationary flat Friedmann universe with dark energy, 
which may support the claim in ref. This can be ex- 
plained by the fact that the pressure gradient of dark 
energy is in the opposite direction to the density gradi- 
ent. It therefore gives an attractive force which pulls the 
matter inwards, whereas the gravitational force is rcpul- 



t 




FIG. 18: The t-r plane of the cosmological white hole solu- 
tion in the case without a similarity horizon, r = 0'^ and 
-|-oo correspond to the Friedmann and quasi-static infinities, 
respectively. There is a wormhole throat a.t z = z% in the 
positive z region, z = z* (< 0) with < t < oo corresponds 
to a spacelike singularity with positive mass, represented by 
a zigzag line. 



z — -|-0, t — oo 




FIG. 19: The fictitious conformal diagram of the Friedmann- 
constant-velocity wormhole solution. There is a similarity 
horizon at 2: = zi (> 0), corresponding to the cosmological 
event horizon, z = 00 with < t < 00 corresponds to a 
timelike boundary of spacetime at infinity. 



sive. So if the black hole is small and the density gradient 
steep, it can accrete the surrounding mass effectively. 

For 7 = 1/3, we have found numerically that the ratio 
of the size of the black hole event horizon to the Hubble 
length is between and 0.70. This means that the size of 
such a self-similar cosmological black hole has an upper 
limit but it can be arbitrarily small, so black holes which 
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TABLE I: Self-similar solutions with 7 = 1/3 which are 
asymptotic to a Friedmann universe at large distance. The 
numerical values are given by ai ~ —0.0253, a2 — —0.0404, 
Qfa ~ —0.0780 and 04 ~ —0.108 for the choice of the gauge 
constants ao = ^'o = 1- The terms similarity horizon, 
trapping horizon, Friedmann, quasi-static, quasi-Kantowski- 
Sachs, positive-mass singular, negative-mass singular are ab- 
breviated as SH, TH, F, QS, QKS, PMS and NMS, respec- 
tively. The last two columns give the number of similarity 
and trapping horizons, with 1(d) indicating one degenerate 
horizon. We note that the spacetime for Aa — ag has not 
been determined. 





Spacetime 


Asymptote 


# SH 


# TH 




Naked singularity 


F-NMS 


1 


1 





Friedmann universe 


F 


1 


1 




Black hole 


F-QKS-PMS 


2 


2 


Qfl 


Black hole 


F-QKS-PMS 


2 


1 (d) 




Wormholo 


F-QKS-QF 


2 





Q!2 


Wormhole 


F-QKS-F 


2 







Wormhole 


F-QKS-QF 


2 





03 


? 


7 


? 


7 




White hole 


F-QS-PMS 


2 





Cti 


White hole 


F-QS-PMS 


1(d) 







White hole 


F-QS-PMS 









These solutions provide intriguing examples of dynam- 
ical wormlioles. Since they are cosmological solutions, 
they are not asymptotically flat and they do not have 
a regular past null infinity. Rather they are asymptotic 
to models which are expanding at spacelike infinity and 
have an initial singularity. 

In such situations a wormhole throat should be locally 
defined as a two-sphere of minimal area on the spacelike 
hypersurfaces rather than a trapping horizon. According 
to the conventional definition of traversability, none of 
our numerical solutions represents a traversable worm- 
hole because of the absence of a past null infinity. The 
strong energy condition is violated there but the null, 
dominant and weak energy conditions are still all sat- 
isfied. Our wormhole solutions do not have a trapping 
horizon: because there is no violation of the null energy 
condition, they are not dynamical wormholes in the sense 
of Hochberg and Visser [l^ or Hayward [3 • 

Finally, we note that the homothetic (rather than co- 
moving) approach is another powerful method of study- 
ing self-similar solutions. In the positive pressure case, 
the homothetic and comoving approaches are comple- 
mentary, which leads to a comprehensive understanding 
of self-similar solutions [l^, [13, HI] ■ A dynamical systems 
approach in the negative pressure case should also shed 
light on various aspects of the problem. 



grow as fast as the universe cannot be too large. It should 
be stressed that there is a one-parameter family of self- 
similar cosmological black hole solutions, so we do not 
have to fine-tune the mass of the black hole to get self- 
similar growth. This means that self-similar black hole 
growth is very plausible in a universe with dark energy 
of this kind. This is in contrast to the usual Newtonian 
argument for a positive pressure fiuid 0, H^l • 

In addition to the black hole solutions, we have found 
self-similar solutions which represent a Friedmann uni- 
verse connected to cither another Friedmann universe or 
some other cosmological model. They are interpreted 
as self-similar cosmological white hole or wormhole so- 
lutions. These are classified into three types, accord- 
ing to whether the other infinity is Friedmann, quasi- 
Friedmann, or quasi-static, among which the last type is 
interpreted as a white hole in a Friedmann background. 
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